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MAGNETIC  TURBULENCE  IN  CLUSTERS  OF  GALAXIES 

T.  A.  EnBlin,1  T.  Clarke,2  C.  Vogt,3  A.  Waelkens,1  A.  A.  Schekochihin4 

RESUMEN 

Los  cumulos  de  galaxias  son  grandes  laboratories  de  turbulencia  en  plasmas  que  nos  permiten  confrontar 
conceptos  teoricos  del  origen  del  campo  magnetico  con  observaciones  detalladas.  La  turbulencia  magnetica  en 
estos  cumulos  puede  ser  estudiada  a  traves  de  la  emision  radio-sincrotron  del  medio  intra-ciimulo  en  forma  de 
relicas  y  halos  de  los  cumulos.  El  espectro  de  potencias  de  los  campos  magneticos  turbulentos  se  puede  estudiar 
a  traves  del  analisis  de  la  rotation  de  Faraday  de  las  fuentes  de  radio  extendidas.  En  el  caso  del  nucleo  frio 
Hydra  A,  el  espectro  magnetico  observado  se  puede  entender  en  terminos  de  una  retroalimentacion,  mediada 
por  la  turbulencia,  entre  el  enfriamiento  del  gas  y  la  actividad  del  chorro  de  la  galaxia  central.  Finalmente 
se  discuten  metodos  para  medir  estadfsticas  de  ordenes  superiores  del  campo  magnetico  usando  correlaciones 
de  parametros  de  Stokes,  los  cuales  nos  permiten  determinar  el  espectro  de  potencias  de  la  fuerza  de  tension 
magnetica.  Dicha  cantidad  estadfstica  de  cuarto  orden  ofrece  una  forma  para  discriminar  diferentes  escenarios 
de  turbulencia  magnetica,  y  diferentes  estructuras  del  campo  magnetico  mediante  observaciones  polarimetricas 
en  radio. 


ABSTRACT 

Galaxy  clusters  are  large  laboratories  for  magnetic  plasma  turbulence  which  permit  us  to  confront  our  theo¬ 
retical  concepts  of  magnetogenesis  with  detailed  observations.  Magnetic  turbulence  in  clusters  can  be  studied 
via  the  radio-synchrotron  emission  from  the  intra-cluster  medium  in  the  form  of  cluster  radio  relics  and  halos. 
The  power  spectrum  of  turbulent  magnetic  fields  can  be  examined  via  Faraday  rotation  analysis  of  extended 
radio  sources.  In  case  of  the  Hydra  A  cool  core,  the  observed  magnetic  spectrum  can  be  understood  in  terms 
of  a  turbulence-mediated  feedback  loop  between  gas  cooling  and  the  jet  activity  of  the  central  galaxy.  Finally, 
methods  to  measure  higher-order  statistics  of  the  magnetic  field  using  Stokes-parameter  correlations  are  dis¬ 
cussed,  which  permit  us  to  determine  the  power  spectrum  of  the  magnetic  tension  force.  This  fourth-order 
statistical  quantity  offers  a  way  to  discriminate  between  different  magnetic  turbulence  scenarios  and  different 
field  structures  using  radio  polarimetric  observations. 

Key  Words:  galaxies:  clusters:  general  -  -  magnetic  fields  —  turbulence 


1.  INTRODUCTION 

Clusters  of  galaxies  are  the  largest  known  mag¬ 
netized  objects  in  the  Universe.  They  are  excellent 
laboratories  to  study  the  interaction  and  symbiosis 
of  turbulent  fluid  motions  and  magnetic  fields.  Mag¬ 
netic  fields  in  clusters  can  be  observed  via  the  radio 
synchrotron  emission  of  relativistic  electron  popula¬ 
tions,  which  can  be  found  predominately  in  merg¬ 
ing  clusters.  Radio  polarimetry  opens  the  possibil¬ 
ity  to  study  the  magnetic  field  properties  in  great 
detail,  since  the  emission’s  intrinsic  polarisation  con¬ 
tains  information  of  the  magnetic  field  orientation  in 
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the  plane  of  the  sky,  and  the  Faraday  rotation  im¬ 
printed  by  the  magnetized  plasma  of  the  intra-cluster 
medium  (ICM)  reveals  the  magnetic  field  component 
along  the  line  of  sight.  Unfortunately,  no  full  3-D  ob¬ 
servable  of  magnetic  fields  is  available  today. 

The  question  of  how  the  theoretical  concepts  of 
ICM  magnetogenesis  and  the  observational  data  can 
be  compared  arises.  In  most  cases  this  has  to  hap¬ 
pen  in  data  space,  since  the  observational  process 
erases  much  of  the  information  of  the  magnetic  fields, 
whereas  synthetic  observations  of  simulated  fields 
should  always  be  possible.  In  a  few  cases  the  compar¬ 
ison  can  be  done  directly  in  theory  space,  in  terms  of 
magnetic  properties,  like  the  power  spectrum.  Usu¬ 
ally,  observations  of  spatially  extended  emission  re¬ 
gions  in  conjunction  with  more  or  less  plausible  as¬ 
sumptions  to  fill  in  missing  information  are  necessary 
for  this  to  be  possible.  Since  such  observational  re¬ 
sults  rely  and  depend  on  the  assumptions  made  (sta- 
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Fig.  1.  Diffuse  radio  emission  of  Abell  2256.  A  Mpc-sized  radio  relic  marks  the  probable  location  of  a  cluster  merger 
shock  wave  (top-right  structure),  and  a  radio  halo  indicates  the  existence  of  turbulent  magnetic  fields  in  the  cluster 
center  (roundish  central  object).  From  Clarke  &  Enfilin  (2006). 


tistical  isotropy,  spherical  symmetry  of  ICM  density, 
etc.),  any  physical  interpretation  has  to  be  aware  of 
this. 

The  next  question  is  in  which  theoretical  frame¬ 
work  the  observational  data  should  be  interpreted. 
Here,  a  simplistic  scenario  for  magnetogenesis  is 
adapted  and  applied  to  observations.  Namely,  it  is 
assumed  that  the  observed  magnetic  fields  are  gen¬ 
erated  and  maintained  by  a  small-scale  dynamo.  It 
is  driven  by  the  turbulence,  which,  initially  hydro- 
dynamic,  becomes  magnetohydrodynamic  when  the 
small-scale  dynamo  saturates.  In  this  picture,  the 
magnetic  fields  strength  and  length-scale  are  usu¬ 
ally  related  to  the  level  of  kinetic  energy  and  turbu¬ 
lence  injection  scale,  respectively.  In  the  description 
of  Subramanian,  Shukurov,  &  Haugen  (2006)  of  the 


saturated  state  of  such  a  dynamo,  the  magnetic  en¬ 
ergy  density  and  length  scales  are  both  smaller  than 
their  hydrodynamical  counterparts,  by  some  factors 
related  to  the  critical  magnetic  Reynolds  number  for 
small-scale  dynamo.  Note  that  the  alternative  pic¬ 
ture  for  the  saturation  of  the  small-scale  dynamo  de¬ 
veloped  by  Schekochihin  et  al.  (2004)  on  the  basis 
of  their  MHD  numerical  simulations  envisions  com¬ 
parable  magnetic  and  kinetic  energies  and  magnetic 
fields  organized  in  long  folded  sheets  with  field  rever¬ 
sal  scale  determined  by  the  dissipation  physics  (the 
true  magnetic  Reynolds  number).  However,  since 
the  MHD  description  fails  at  small  scales  in  weakly 
collisional  cluster  plasmas,  one  has  to  invoke  some 
yet  poorly  understood  plasma  micro-physics  mecha¬ 
nisms  to  determine  the  effective  magnetic  Reynolds 
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number  for  the  field  reversals  (Schekochihin  &  Cow¬ 
ley  2006).  Here  we  do  not  discuss  in  detail  these  var¬ 
ious  scenarios  and  the  uncertainties  associated  with 
them,  but  instead  adopt  the  view  that  the  magnetic 
turbulence  in  clusters  can  be  parametrized  by  some 
effective  magnetic  Reynolds  number  of  order  ~  102 
(see  further  details  in  Enfilin  et  al.  2006).  The  bot¬ 
tom  line  of  this  talk  is  that  this  simplistic  picture 
seems  to  be  in  good  agreement  with  the  observations 
so  far.  Novel  methods  to  extract  more  information 
on  fields  from  observations,  in  order  to  test  this  and 
other  pictures  further,  are  currently  under  develop¬ 
ment  are  are  presented  at  the  end. 

2.  CLUSTER  RADIO  HALOS  AND  RELICS 

Many  merging  clusters  of  galaxies  exhibit  Mpc- 
sized  synchrotron  radio  emission  regions.  These  can 
be  classified  into  two  categories: 

-  Cluster  radio  halos  are  cluster  wide,  unpolar¬ 
ized  emission  regions.  Their  morphology  strongly  re¬ 
sembles  the  X-ray  emitting  gas  morphology  of  clus¬ 
ters,  indicating  that  they  are  volume-probes  of  mag¬ 
netic  fields  and  relativistic  electrons.  The  origin  of 
the  electrons  is  unclear,  they  could  be  either  main¬ 
tained  by  Fermi-II  re-acceleration  (Giovannini  et  al. 
1993;  Brunetti  et  al.  2001),  or  freshly  injected  by 
hadronic  interactions  of  a  long-lived  relativistic  pro¬ 
ton  population  (Dennison  1980). 

-  Cluster  radio  relics  are  typically  peripheral,  po¬ 
larized  emission  regions,  which  sometimes  exhibit  fil¬ 
amentary  structures.  At  least  the  Mpc-sized  relics 
are  probalbly  due  to  Fermi-I  acceleration  of  electrons 
at  accretion  or  merger  shock  waves  (Enfilin  et  al. 
1998;  Giacintucci  et  al.  2008).  Since  the  electrons 
cool  rapidly  within  only  10s  years,  they  illuminate 
the  magnetic  fields  only  in  sheet-like  sub- volumes  of 
the  clusters. 

An  example  of  a  cluster  showing  both  phenom¬ 
ena,  relics  and  halos,  is  given  in  Figure  1.  There, 
Mpc-sized,  filamentary  magnetic  structures  can  be 
seen  within  the  relic  region,  very  similar  to  the  sheet 
like  structures  seen  in  slices  through  magnetic  tur¬ 
bulence  simulations  in  Schekochihin  et  al.  (2004),  see 
also  Figure  3. 

It  is  quite  obvious  why  only  merging  cluster  ex¬ 
hibit  radio  relics,  since  strong  shock  waves  can  only 
be  found  there.  The  rareness  of  radio  halos  lies  cer¬ 
tainly  partly  in  the  fact,  that  only  clusters  with  well 
developed  turbulence  may  be  able  to  maintain  suf¬ 
ficiently  strong  magnetic  fields  by  a  small  scale  dy¬ 
namo. 

In  any  case,  cluster  wide  radio  emission  regions 
like  halos  and  relics  permit  us  to  probe  magnetic 


Fig.  2.  The  magnetic  power  spectra  of  the  cool  core  of 
the  Hydra  cluster  Vogt  &  Enfilin  (2005). 

fields  on  the  largest  scales  accessible  so  far.  Halos  are 
best  suited  to  provide  us  with  a  global  view  on  the 
field  distribution  in  a  galaxy  cluster,  whereas  relics 
highlight  small-scale  details  and  the  geometry  of  in¬ 
dividual  flux  bundles. 

3.  MAGNETIC  TURBULENCE  IN  COOL  CORES 

Relaxed  clusters  also  exhibit  magnetic  fields, 
where  the  strongest  fields  are  actually  found  in  the 
most  relaxed  clusters  which  developed  a  cool  core 
structure.  Cluster  cool  cores  are  central  regions  of 
clusters  in  which  a  cooling  instability  of  the  gas  has 
developed.  The  strong  energy  losses  due  to  X-ray 
emission  of  the  dense  central  gas  are  partly  com¬ 
pensated  by  an  efficient  energy  feedback  process.  It 
is  believed  that  gas  condensing  from  the  ICM  and 
falling  into  the  central  galaxy  feeds  the  central  black 
hole.  The  black  hole  becomes  thereby  active  and 
emits  powerful  radio  jets.  The  radio  bubbles  inflated 
by  these  jets  drive  weak  shocks  into  the  cool  core  gas, 
and  stir  turbulence  during  their  buoyant  movement 
in  the  cluster  atmosphere.  This  turbulence,  on  the 
one  hand,  heats  the  medium,  and  counteracts  the 
cooling  instability,  thus  providing  a  negative  feed¬ 
back.  On  the  other  hand,  this  turbulence  should 
maintain  some  level  of  magnetization,  and  thereby 
imprint  its  characteristic  properties  (energy  density, 
correlation  length)  onto  the  magnetic  field. 

In  the  simplistic  picture  of  magnetogenesis 
adopted  here,  the  turbulence  maintains  magnetic 
fields  on  length  scales  somewhat  smaller  than  its  own 
driving  scale,  which  is  the  size  of  the  radio  bubbles. 
Thus  field  correlation  length  of  a  few  kpc,  and  field 
strength  of  the  order  of  ~  5  are  expected  in  a 
typical  cool  core.  The  ICM  fields  seem  also  to  play 
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Fig.  3.  Magnetic  energy  density  in  a  MHD  simulation  (left,  the  data  was  taken  from  run  S5  of  Schekochihin  et  al.  2004) 
and  in  a  Gaussian  random  field  setup  (right). 
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Fig.  4.  Magnetic  tension  force  density  in  a  MHD  simulation  (left,  the  data  was  taken  from  run  S5  of  Schekochihin  et 
al.  2004)  and  in  a  Gaussian  random  field  setup  (right). 
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some  role  in  stabilizing  the  rising  bubbles  against 
disruptive  instabilities  (Ruszkowski  et  al.  2008). 

A  Faraday  rotation  based  determination  of  the 
magnetic  field  strength  could  actually  confirm  these 
expectations.  In  the  case  of  the  Hydra  A  cluster, 
Vogt  &  EnBlin  (2005)  measured  the  magnetic  power 
spectrum,  as  shown  in  Figure  2,  while  assuming 
statistically  isotropic  magnetic  fields  and  that  the 
mean  magnetic  energy  density  scales  linearly  with 


that  of  the  cluster  gas  atmosphere.  The  observed 
Kolmogorov-like  power  law  indicates  a  turbulent  ori¬ 
gin  of  the  fields.  The  field  strength  of  7  /rG  and 
length  scale  of  3  kpc  fit  well  with  the  expectations 
for  a  steady  state  turbulence  level  mediating  the  en¬ 
ergy  feedback  between  the  radio  bubbles  and  the  X- 
ray  cooling  cool  core  gas  (EnBlin  &  Vogt  2006).  A 
similar  consistency  between  the  Faraday  signature 
of  cool  core  magnetic  fields  and  their  strength  and 
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Fig.  5.  Left:  The  magnetic  power  spectra  of  the  two  field  configurations  are  very  similar,  although  the  visual  inspection  of 
the  fields  reveals  significant  differences.  Right:  The  magnetic  tension  force  power  spectrum  of  the  two  field  configurations 
differ  substantially,  and  thereby  permit  to  discriminate  the  different  turbulence  scenarios. 


length  scale  expected  in  this  scenario  could  also  be 
demonstrated  for  a  number  of  other  clusters  (Enfilin 
&  Vogt  2006). 

4.  MEASURING  THE  MAGNETIC  TENSION 
FORCE  SPECTRUM 


A  detailed  investigation  of  magnetic  field  statis¬ 
tics  in  galaxy  clusters  should  also  try  to  extract 
higher-order  correlation  functions.  Specially  mean¬ 
ingful  would  be  a  measurement  of  the  tension  force 
power  spectrum.  The  tension  force  is  F  =  B  •  VB, 
where  B  is  the  magnetic  field.  This  quantity  mea¬ 
sures  how  strong  the  magnetic  back-reaction  forces 
are  localized  or  distributed  through  out  the  volume. 
It  also  conveys  information  about  the  geometrical 
structure  of  the  magnetic  fields  because  it  involves 
the  field  gradient  along  the  field’s  own  local  direc¬ 
tion.  In  particular,  measuring  the  tension  force 
spectrum  allows  one  to  distinguish  folded  filamented 
magnetic  fields  (such  as  found  in  numerical  simu¬ 
lations)  from  isotropically  tangled  fields  sometimes 
naively  assumed  on  the  basis  of  the  power  spectra 
alone.  Thus,  statistics  of  the  tension  force  are  very 
discriminative  of  magnetic  turbulence  scenarios. 

This  is  illustrated  by  Figure  3,  where  a  field  con¬ 
figuration  of  a  MHD-turbulence  simulation  is  con¬ 
fronted  to  divergence-free  Gaussian  random  fields 
with  the  same  magnetic  power  spectra,  as  shown 
in  Figure  5.  The  spatial  distribution  of  the  tension 
forces  in  the  two  scenarios  is  presented  in  Figure  4. 
The  two  tension  force  spectra  differ,  as  Figure  5  con¬ 
firms,  in  that  the  forces  are  stronger  and  more  con¬ 
centrated  on  small  scales  in  the  Gaussian  random 
held  as  compared  to  the  simulation. 

If  one  assumes  that  the  magnetic  fields  are  radio¬ 
synchrotron  illuminated  by  a  spatially  homogeneous 


cosmic  ray  electron  population,  observables  on  the 
magnetic  fields  are  available.  The  polarized  radio 
intensity  emitted  by  the  synchrotron  process  is,  sim¬ 
ilarly  to  the  tension  force,  quadratically  dependent 
on  the  magnetic  held  strength.  Thus,  its  spectrum  is 
a  fourth  order  function  of  the  magnetic  held,  as  the 
tension  force  spectrum  is.  What  is  not  clear  ab  initio 
is  if  the  polarization  signal  contains  all  information 
in  order  to  reconstruct  the  tension  spectrum.  Actu¬ 
ally,  it  can  be  shown  that  this  is  not  the  case,  since 
the  synchrotron  emission  is  insensitive  to  the  line-of 
sight  component  of  the  held  which  enters,  however, 
in  the  tension  spectrum.  This  missing  information 
can  only  be  replaced  by  a  plausible  assumption,  e.g. 
in  our  work  the  assumption  of  statistical  isotropy  of 
the  helds. 

Polarised  emissivity  maps  are  shown  in  Figure  6. 
Waelkens  et  al.  (in  prep.)  -see  also  Enhlin  et  al. 
(2006)-  managed  to  demonstrate  that  the  cross-  and 
auto-  correlation  functions  of  the  Stokes  I,  Q  and 
U  parameters  of  the  synchrotron  emission  contain 
enough  information  to  reconstruct  the  tension  force 
power  spectrum.  This  is  possible  under  the  assump¬ 
tions  of  statistically  homogeneous  and  isotropic  mag¬ 
netic  helds,  and  a  spatially  homogeneous  cosmic  ray 
electron  distribution  with  a  spectral  index  close  to 
—3.  In  such  circumstances,  which  are  not  too  arti¬ 
ficial  in  many  environments,  the  tension  force  spec¬ 
trum  can  be  measured,  and  the  measured  spectra 
indeed  permit  us  to  discriminate  the  two  magnetic 
turbulence  scenarios.  This  is  illustrated  by  Figure  7. 

Since  Faraday  rotation  suppresses  the  radio  halo 
polarisation  signal,  cluster  radio  relics  are  probably 
a  likely  place  where  this  method  can  be  applied.  It 
should  also  be  possible  to  investigate  the  magnetic 
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Fig.  6.  Polarized  intensity  of  the  radio  synchrotron  emission  from  the  magnetic  fields  of  the  MHD  simulation  (left,  the 
data  was  taken  from  run  S5  of  Schekochihin  et  al.  2004)  and  of  the  Gaussian  random  field  setup  (right). 


Fig.  7.  Magnetic  tension  force  spectrum  recovered  from  the  polarized  intensity  in  comparison  to  the  correct  one.  The 
case  of  magnetic  fields  of  the  MHD  simulation  is  shown  on  the  left  and  that  of  the  Gaussian  random  field  setup  on  the 
right. 


fields  inside  the  radio  bubbles,  since  these  seem  to 
be  Faraday  rotation  free. 


cluster  medium  should  become  a  very  good  test  bed 
for  theories  of  turbulent  magnetogenesis. 


5.  OUTLOOK 

With  the  upcoming  sensitive  radio  telescopes  of 
the  next  generation  like  LOFAR,  LWA,  eVLA,  and 
SKA,  the  number  of  detected  cluster  radio  halos  and 
relics  should  strongly  increase  (EnBlin  &  Rottgering 
2002;  Briiggen  et  al.  2003;  Cassano  &  Brunetti  2005). 
Also  many  more  details  of  the  fields  should  be¬ 
come  visible  due  to  higher  sensitivity,  resolution  and 
the  development  of  new  observational  techniques, 
like  the  Faraday  tomography.  Therefore,  the  intra¬ 


TE  thanks  the  SOC  for  the  invitation  to  this  ex¬ 
cellent  conference.  AW  thanks  Francisco  Kitaura  for 
discussion  on  power  spectrum  realization. 
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